home *** CD-ROM | disk | FTP | other *** search
Text File | 2003-02-09 | 60.8 KB | 1,072 lines |
- ;;; Compiled by f2cl version 2.0 beta 2002-05-06
- ;;;
- ;;; Options: ((:prune-labels nil) (:auto-save t) (:relaxed-array-decls t)
- ;;; (:coerce-assigns :as-needed) (:array-type ':simple-array)
- ;;; (:array-slicing nil) (:declare-common nil)
- ;;; (:float-format double-float))
-
- (in-package "SLATEC")
-
-
- (let ((ar (make-array 14 :element-type 'double-float))
- (br (make-array 14 :element-type 'double-float))
- (c (make-array 105 :element-type 'double-float))
- (alfa (make-array 180 :element-type 'double-float))
- (beta (make-array 210 :element-type 'double-float))
- (gama (make-array 30 :element-type 'double-float))
- (ex1 0.33333333333333337)
- (ex2 0.6666666666666667)
- (hpi 1.5707963267948966)
- (gpi 3.141592653589793)
- (thpi 4.71238898038469)
- (zeror 0.0)
- (zeroi 0.0)
- (coner 1.0)
- (conei 0.0))
- (declare (type double-float conei coner zeroi zeror thpi gpi hpi ex2 ex1)
- (type (simple-array double-float (30)) gama)
- (type (simple-array double-float (210)) beta)
- (type (simple-array double-float (180)) alfa)
- (type (simple-array double-float (105)) c)
- (type (simple-array double-float (14)) br ar))
- (f2cl-lib:fset (f2cl-lib:fref ar (1) ((1 14))) 1.0)
- (f2cl-lib:fset (f2cl-lib:fref ar (2) ((1 14))) 0.10416666666666669)
- (f2cl-lib:fset (f2cl-lib:fref ar (3) ((1 14))) 0.08355034722222221)
- (f2cl-lib:fset (f2cl-lib:fref ar (4) ((1 14))) 0.12822657455632716)
- (f2cl-lib:fset (f2cl-lib:fref ar (5) ((1 14))) 0.29184902646414046)
- (f2cl-lib:fset (f2cl-lib:fref ar (6) ((1 14))) 0.8816272674437576)
- (f2cl-lib:fset (f2cl-lib:fref ar (7) ((1 14))) 3.3214082818627677)
- (f2cl-lib:fset (f2cl-lib:fref ar (8) ((1 14))) 14.995762986862555)
- (f2cl-lib:fset (f2cl-lib:fref ar (9) ((1 14))) 78.92301301158652)
- (f2cl-lib:fset (f2cl-lib:fref ar (10) ((1 14))) 474.4515388682643)
- (f2cl-lib:fset (f2cl-lib:fref ar (11) ((1 14))) 3207.490090890662)
- (f2cl-lib:fset (f2cl-lib:fref ar (12) ((1 14))) 24086.549640874004)
- (f2cl-lib:fset (f2cl-lib:fref ar (13) ((1 14))) 198923.1191695098)
- (f2cl-lib:fset (f2cl-lib:fref ar (14) ((1 14))) 1791902.0077753437)
- (f2cl-lib:fset (f2cl-lib:fref br (1) ((1 14))) 1.0)
- (f2cl-lib:fset (f2cl-lib:fref br (2) ((1 14))) -0.14583333333333334)
- (f2cl-lib:fset (f2cl-lib:fref br (3) ((1 14))) -0.09874131944444445)
- (f2cl-lib:fset (f2cl-lib:fref br (4) ((1 14))) -0.14331205391589505)
- (f2cl-lib:fset (f2cl-lib:fref br (5) ((1 14))) -0.31722720267841353)
- (f2cl-lib:fset (f2cl-lib:fref br (6) ((1 14))) -0.9424291479571203)
- (f2cl-lib:fset (f2cl-lib:fref br (7) ((1 14))) -3.5112030408263544)
- (f2cl-lib:fset (f2cl-lib:fref br (8) ((1 14))) -15.727263620368046)
- (f2cl-lib:fset (f2cl-lib:fref br (9) ((1 14))) -82.28143909718595)
- (f2cl-lib:fset (f2cl-lib:fref br (10) ((1 14))) -492.35537052367056)
- (f2cl-lib:fset (f2cl-lib:fref br (11) ((1 14))) -3316.218568547972)
- (f2cl-lib:fset (f2cl-lib:fref br (12) ((1 14))) -24827.67424520859)
- (f2cl-lib:fset (f2cl-lib:fref br (13) ((1 14))) -204526.58731512981)
- (f2cl-lib:fset (f2cl-lib:fref br (14) ((1 14))) -1838444.91706821)
- (f2cl-lib:fset (f2cl-lib:fref c (1) ((1 105))) 1.0)
- (f2cl-lib:fset (f2cl-lib:fref c (2) ((1 105))) -0.20833333333333337)
- (f2cl-lib:fset (f2cl-lib:fref c (3) ((1 105))) 0.125)
- (f2cl-lib:fset (f2cl-lib:fref c (4) ((1 105))) 0.3342013888888889)
- (f2cl-lib:fset (f2cl-lib:fref c (5) ((1 105))) -0.40104166666666674)
- (f2cl-lib:fset (f2cl-lib:fref c (6) ((1 105))) 0.0703125)
- (f2cl-lib:fset (f2cl-lib:fref c (7) ((1 105))) -1.0258125964506173)
- (f2cl-lib:fset (f2cl-lib:fref c (8) ((1 105))) 1.8464626736111112)
- (f2cl-lib:fset (f2cl-lib:fref c (9) ((1 105))) -0.8912109375)
- (f2cl-lib:fset (f2cl-lib:fref c (10) ((1 105))) 0.0732421875)
- (f2cl-lib:fset (f2cl-lib:fref c (11) ((1 105))) 4.669584423426247)
- (f2cl-lib:fset (f2cl-lib:fref c (12) ((1 105))) -11.207002616222994)
- (f2cl-lib:fset (f2cl-lib:fref c (13) ((1 105))) 8.78912353515625)
- (f2cl-lib:fset (f2cl-lib:fref c (14) ((1 105))) -2.3640869140625)
- (f2cl-lib:fset (f2cl-lib:fref c (15) ((1 105))) 0.112152099609375)
- (f2cl-lib:fset (f2cl-lib:fref c (16) ((1 105))) -28.212072558200244)
- (f2cl-lib:fset (f2cl-lib:fref c (17) ((1 105))) 84.63621767460074)
- (f2cl-lib:fset (f2cl-lib:fref c (18) ((1 105))) -91.81824154324002)
- (f2cl-lib:fset (f2cl-lib:fref c (19) ((1 105))) 42.53499874538846)
- (f2cl-lib:fset (f2cl-lib:fref c (20) ((1 105))) -7.368794359479632)
- (f2cl-lib:fset (f2cl-lib:fref c (21) ((1 105))) 0.22710800170898437)
- (f2cl-lib:fset (f2cl-lib:fref c (22) ((1 105))) 212.57013003921713)
- (f2cl-lib:fset (f2cl-lib:fref c (23) ((1 105))) -765.2524681411817)
- (f2cl-lib:fset (f2cl-lib:fref c (24) ((1 105))) 1059.9904525279999)
- (f2cl-lib:fset (f2cl-lib:fref c (25) ((1 105))) -699.5796273761325)
- (f2cl-lib:fset (f2cl-lib:fref c (26) ((1 105))) 218.1905117442116)
- (f2cl-lib:fset (f2cl-lib:fref c (27) ((1 105))) -26.491430486951554)
- (f2cl-lib:fset (f2cl-lib:fref c (28) ((1 105))) 0.5725014209747314)
- (f2cl-lib:fset (f2cl-lib:fref c (29) ((1 105))) -1919.457662318407)
- (f2cl-lib:fset (f2cl-lib:fref c (30) ((1 105))) 8061.722181737309)
- (f2cl-lib:fset (f2cl-lib:fref c (31) ((1 105))) -13586.550006434138)
- (f2cl-lib:fset (f2cl-lib:fref c (32) ((1 105))) 11655.393336864532)
- (f2cl-lib:fset (f2cl-lib:fref c (33) ((1 105))) -5305.646978613403)
- (f2cl-lib:fset (f2cl-lib:fref c (34) ((1 105))) 1200.9029132163525)
- (f2cl-lib:fset (f2cl-lib:fref c (35) ((1 105))) -108.09091978839464)
- (f2cl-lib:fset (f2cl-lib:fref c (36) ((1 105))) 1.7277275025844574)
- (f2cl-lib:fset (f2cl-lib:fref c (37) ((1 105))) 20204.291330966145)
- (f2cl-lib:fset (f2cl-lib:fref c (38) ((1 105))) -96980.5983886375)
- (f2cl-lib:fset (f2cl-lib:fref c (39) ((1 105))) 192547.00123253153)
- (f2cl-lib:fset (f2cl-lib:fref c (40) ((1 105))) -203400.17728041552)
- (f2cl-lib:fset (f2cl-lib:fref c (41) ((1 105))) 122200.46498301746)
- (f2cl-lib:fset (f2cl-lib:fref c (42) ((1 105))) -41192.65496889756)
- (f2cl-lib:fset (f2cl-lib:fref c (43) ((1 105))) 7109.514302489364)
- (f2cl-lib:fset (f2cl-lib:fref c (44) ((1 105))) -493.915304773088)
- (f2cl-lib:fset (f2cl-lib:fref c (45) ((1 105))) 6.074042001273483)
- (f2cl-lib:fset (f2cl-lib:fref c (46) ((1 105))) -242919.1879005513)
- (f2cl-lib:fset (f2cl-lib:fref c (47) ((1 105))) 1311763.6146629772)
- (f2cl-lib:fset (f2cl-lib:fref c (48) ((1 105))) -2998015.9185381066)
- (f2cl-lib:fset (f2cl-lib:fref c (49) ((1 105))) 3763271.297656404)
- (f2cl-lib:fset (f2cl-lib:fref c (50) ((1 105))) -2813563.226586534)
- (f2cl-lib:fset (f2cl-lib:fref c (51) ((1 105))) 1268365.2733216248)
- (f2cl-lib:fset (f2cl-lib:fref c (52) ((1 105))) -331645.1724845636)
- (f2cl-lib:fset (f2cl-lib:fref c (53) ((1 105))) 45218.76898136273)
- (f2cl-lib:fset (f2cl-lib:fref c (54) ((1 105))) -2499.8304818112097)
- (f2cl-lib:fset (f2cl-lib:fref c (55) ((1 105))) 24.380529699556064)
- (f2cl-lib:fset (f2cl-lib:fref c (56) ((1 105))) 3284469.853072038)
- (f2cl-lib:fset (f2cl-lib:fref c (57) ((1 105))) -1.9706819118432228e+7)
- (f2cl-lib:fset (f2cl-lib:fref c (58) ((1 105))) 5.095260249266464e+7)
- (f2cl-lib:fset (f2cl-lib:fref c (59) ((1 105))) -7.410514821153266e+7)
- (f2cl-lib:fset (f2cl-lib:fref c (60) ((1 105))) 6.634451227472903e+7)
- (f2cl-lib:fset (f2cl-lib:fref c (61) ((1 105))) -3.756717666076335e+7)
- (f2cl-lib:fset (f2cl-lib:fref c (62) ((1 105))) 1.3288767166421817e+7)
- (f2cl-lib:fset (f2cl-lib:fref c (63) ((1 105))) -2785618.1280864547)
- (f2cl-lib:fset (f2cl-lib:fref c (64) ((1 105))) 308186.4046126624)
- (f2cl-lib:fset (f2cl-lib:fref c (65) ((1 105))) -13886.089753717042)
- (f2cl-lib:fset (f2cl-lib:fref c (66) ((1 105))) 110.01714026924674)
- (f2cl-lib:fset (f2cl-lib:fref c (67) ((1 105))) -4.932925366450996e+7)
- (f2cl-lib:fset (f2cl-lib:fref c (68) ((1 105))) 3.2557307418576575e+8)
- (f2cl-lib:fset (f2cl-lib:fref c (69) ((1 105))) -9.394623596815784e+8)
- (f2cl-lib:fset (f2cl-lib:fref c (70) ((1 105))) 1.55359689957058e+9)
- (f2cl-lib:fset (f2cl-lib:fref c (71) ((1 105))) -1.6210805521083374e+9)
- (f2cl-lib:fset (f2cl-lib:fref c (72) ((1 105))) 1.1068428168230145e+9)
- (f2cl-lib:fset (f2cl-lib:fref c (73) ((1 105))) -4.958897842750303e+8)
- (f2cl-lib:fset (f2cl-lib:fref c (74) ((1 105))) 1.4206290779753308e+8)
- (f2cl-lib:fset (f2cl-lib:fref c (75) ((1 105))) -2.4474062725738724e+7)
- (f2cl-lib:fset (f2cl-lib:fref c (76) ((1 105))) 2243768.1779224495)
- (f2cl-lib:fset (f2cl-lib:fref c (77) ((1 105))) -84005.4336030241)
- (f2cl-lib:fset (f2cl-lib:fref c (78) ((1 105))) 551.3358961220206)
- (f2cl-lib:fset (f2cl-lib:fref c (79) ((1 105))) 8.147890961183122e+8)
- (f2cl-lib:fset (f2cl-lib:fref c (80) ((1 105))) -5.866481492051847e+9)
- (f2cl-lib:fset (f2cl-lib:fref c (81) ((1 105))) 1.8688207509295826e+10)
- (f2cl-lib:fset (f2cl-lib:fref c (82) ((1 105))) -3.4632043388158773e+10)
- (f2cl-lib:fset (f2cl-lib:fref c (83) ((1 105))) 4.1280185579753975e+10)
- (f2cl-lib:fset (f2cl-lib:fref c (84) ((1 105))) -3.302659974980072e+10)
- (f2cl-lib:fset (f2cl-lib:fref c (85) ((1 105))) 1.79542137311556e+10)
- (f2cl-lib:fset (f2cl-lib:fref c (86) ((1 105))) -6.563293792619285e+9)
- (f2cl-lib:fset (f2cl-lib:fref c (87) ((1 105))) 1.5592798648792575e+9)
- (f2cl-lib:fset (f2cl-lib:fref c (88) ((1 105))) -2.251056618894153e+8)
- (f2cl-lib:fset (f2cl-lib:fref c (89) ((1 105))) 1.7395107553978165e+7)
- (f2cl-lib:fset (f2cl-lib:fref c (90) ((1 105))) -549842.3275722887)
- (f2cl-lib:fset (f2cl-lib:fref c (91) ((1 105))) 3038.0905109223845)
- (f2cl-lib:fset (f2cl-lib:fref c (92) ((1 105))) -1.4679261247695616e+10)
- (f2cl-lib:fset (f2cl-lib:fref c (93) ((1 105))) 1.1449823773202579e+11)
- (f2cl-lib:fset (f2cl-lib:fref c (94) ((1 105))) -3.990961752244665e+11)
- (f2cl-lib:fset (f2cl-lib:fref c (95) ((1 105))) 8.192186695485775e+11)
- (f2cl-lib:fset (f2cl-lib:fref c (96) ((1 105))) -1.0983751560812234e+12)
- (f2cl-lib:fset (f2cl-lib:fref c (97) ((1 105))) 1.008158106865382e+12)
- (f2cl-lib:fset (f2cl-lib:fref c (98) ((1 105))) -6.453648692453764e+11)
- (f2cl-lib:fset (f2cl-lib:fref c (99) ((1 105))) 2.879006499061506e+11)
- (f2cl-lib:fset (f2cl-lib:fref c (100) ((1 105))) -8.786707217802326e+10)
- (f2cl-lib:fset (f2cl-lib:fref c (101) ((1 105))) 1.763473060683497e+10)
- (f2cl-lib:fset (f2cl-lib:fref c (102) ((1 105))) -2.167164983223795e+9)
- (f2cl-lib:fset (f2cl-lib:fref c (103) ((1 105))) 1.4315787671888897e+8)
- (f2cl-lib:fset (f2cl-lib:fref c (104) ((1 105))) -3871833.4425726123)
- (f2cl-lib:fset (f2cl-lib:fref c (105) ((1 105))) 18257.755474293175)
- (f2cl-lib:fset (f2cl-lib:fref alfa (1) ((1 180))) -0.0044444444444444444)
- (f2cl-lib:fset (f2cl-lib:fref alfa (2) ((1 180))) -9.220779220779222e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (3) ((1 180))) -8.848928848928851e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (4) ((1 180))) 1.6592768783244974e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (5) ((1 180))) 2.4669137274179292e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (6) ((1 180))) 2.6599558934625483e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (7) ((1 180))) 2.6182429706150095e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (8) ((1 180))) 2.4873043734465564e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (9) ((1 180))) 2.327210400832321e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (10) ((1 180))) 2.163624857123651e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (11) ((1 180))) 2.0073885876275238e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (12) ((1 180))) 1.8626763663754517e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (13) ((1 180))) 1.730607759178765e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (14) ((1 180))) 1.6109170592901578e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (15) ((1 180))) 1.5027477416090815e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (16) ((1 180))) 1.405034973912698e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (17) ((1 180))) 1.316688165459228e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (18) ((1 180))) 1.2366744559825327e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (19) ((1 180))) 1.1640527147473789e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (20) ((1 180))) 1.0979829837271338e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (21) ((1 180))) 1.0377241042299283e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (22) ((1 180))) 9.826260783693634e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (23) ((1 180))) 9.321205172495034e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (24) ((1 180))) 8.857108524787118e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (25) ((1 180))) 8.429631057157001e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (26) ((1 180))) 8.034975484077912e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (27) ((1 180))) 7.669813453592074e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (28) ((1 180))) 7.33122157481778e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (29) ((1 180))) 7.016626251631414e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (30) ((1 180))) 6.723756337901602e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (31) ((1 180))) 6.937355413545891e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (32) ((1 180))) 2.3224174518292165e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (33) ((1 180))) -1.4198627355669122e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (34) ((1 180))) -1.1644493167204864e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (35) ((1 180))) -1.508035580530488e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (36) ((1 180))) -1.5512192491809623e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (37) ((1 180))) -1.4680975664646556e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (38) ((1 180))) -1.3381550386749136e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (39) ((1 180))) -1.1974497568425405e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (40) ((1 180))) -1.0618431920797404e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (41) ((1 180))) -9.376995498911946e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (42) ((1 180))) -8.269230455881934e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (43) ((1 180))) -7.293743481552213e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (44) ((1 180))) -6.440423577210163e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (45) ((1 180))) -5.69611566009369e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (46) ((1 180))) -5.047310443035617e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (47) ((1 180))) -4.481348680088828e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (48) ((1 180))) -3.9868872771759895e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (49) ((1 180))) -3.554005329720425e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (50) ((1 180))) -3.174142566090225e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (51) ((1 180))) -2.8399679390417487e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (52) ((1 180))) -2.545227206348706e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (53) ((1 180))) -2.2845929716472457e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (54) ((1 180))) -2.0535275310648063e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (55) ((1 180))) -1.8481621762766612e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (56) ((1 180))) -1.665193300213938e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (57) ((1 180))) -1.501794129801195e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (58) ((1 180))) -1.3555403137904054e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (59) ((1 180))) -1.2243474647385814e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (60) ((1 180))) -1.1064188481130818e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (61) ((1 180))) -3.5421197145774386e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (62) ((1 180))) -1.5616126394515942e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (63) ((1 180))) 3.0446550359493646e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (64) ((1 180))) 1.3019865577324272e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (65) ((1 180))) 1.6747110669971227e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (66) ((1 180))) 1.7022258768359258e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (67) ((1 180))) 1.565014276085947e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (68) ((1 180))) 1.3633917097744513e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (69) ((1 180))) 1.1488669202982513e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (70) ((1 180))) 9.458690930346881e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (71) ((1 180))) 7.644984192508983e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (72) ((1 180))) 6.075703349651974e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (73) ((1 180))) 4.743942992905088e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (74) ((1 180))) 3.627575120053444e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (75) ((1 180))) 2.69939714979225e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (76) ((1 180))) 1.9321093824793925e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (77) ((1 180))) 1.300566747939632e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (78) ((1 180))) 7.826208667444966e-6)
- (f2cl-lib:fset (f2cl-lib:fref alfa (79) ((1 180))) 3.592574858193515e-6)
- (f2cl-lib:fset (f2cl-lib:fref alfa (80) ((1 180))) 1.4404004981425184e-7)
- (f2cl-lib:fset (f2cl-lib:fref alfa (81) ((1 180))) -2.653967696979391e-6)
- (f2cl-lib:fset (f2cl-lib:fref alfa (82) ((1 180))) -4.9134686709848585e-6)
- (f2cl-lib:fset (f2cl-lib:fref alfa (83) ((1 180))) -6.727392960912483e-6)
- (f2cl-lib:fset (f2cl-lib:fref alfa (84) ((1 180))) -8.172693796786579e-6)
- (f2cl-lib:fset (f2cl-lib:fref alfa (85) ((1 180))) -9.313047150935612e-6)
- (f2cl-lib:fset (f2cl-lib:fref alfa (86) ((1 180))) -1.0201141879801645e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (87) ((1 180))) -1.088059625105929e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (88) ((1 180))) -1.1387548150960355e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (89) ((1 180))) -1.1751967567455643e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (90) ((1 180))) -1.1998736487094415e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (91) ((1 180))) 3.781941992017729e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (92) ((1 180))) 2.0247195276181618e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (93) ((1 180))) -6.379385063188625e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (94) ((1 180))) -2.385982306030059e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (95) ((1 180))) -3.109162560273616e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (96) ((1 180))) -3.136801152475764e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (97) ((1 180))) -2.789502737913234e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (98) ((1 180))) -2.2856408261914143e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (99) ((1 180))) -1.7524528034084677e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (100) ((1 180))) -1.2554406306069035e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (101) ((1 180))) -8.229828728202083e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (102) ((1 180))) -4.628607305881165e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (103) ((1 180))) -1.7233430236696232e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (104) ((1 180))) 5.606904823046023e-6)
- (f2cl-lib:fset (f2cl-lib:fref alfa (105) ((1 180))) 2.313954431482868e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (106) ((1 180))) 3.62642745856794e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (107) ((1 180))) 4.5800612449018885e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (108) ((1 180))) 5.245952949591141e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (109) ((1 180))) 5.683962085458154e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (110) ((1 180))) 5.943498203931041e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (111) ((1 180))) 6.064785275784218e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (112) ((1 180))) 6.080239077884365e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (113) ((1 180))) 6.015778945394604e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (114) ((1 180))) 5.891996573446985e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (115) ((1 180))) 5.725158237775931e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (116) ((1 180))) 5.528043755858526e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (117) ((1 180))) 5.310637738028802e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (118) ((1 180))) 5.080693020123257e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (119) ((1 180))) 4.844186476200949e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (120) ((1 180))) 4.6056858160747544e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (121) ((1 180))) -6.911413972882943e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (122) ((1 180))) -4.2997663305887196e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (123) ((1 180))) 1.8306773598003906e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (124) ((1 180))) 6.600881475420142e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (125) ((1 180))) 8.75964969951186e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (126) ((1 180))) 8.773352359582357e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (127) ((1 180))) 7.493695853789907e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (128) ((1 180))) 5.638323297569809e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (129) ((1 180))) 3.680593199714432e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (130) ((1 180))) 1.884645355144556e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (131) ((1 180))) 3.7066305766490415e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (132) ((1 180))) -8.285202202321371e-5)
- (f2cl-lib:fset (f2cl-lib:fref alfa (133) ((1 180))) -1.72751952869173e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (134) ((1 180))) -2.3631487360587297e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (135) ((1 180))) -2.779661506949067e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (136) ((1 180))) -3.020795141554569e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (137) ((1 180))) -3.125947126438201e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (138) ((1 180))) -3.128725587580672e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (139) ((1 180))) -3.056780384663244e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (140) ((1 180))) -2.932264706145573e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (141) ((1 180))) -2.772556555829348e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (142) ((1 180))) -2.5910392846703173e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (143) ((1 180))) -2.3978401439648037e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (144) ((1 180))) -2.2004826004542286e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (145) ((1 180))) -2.0044391109497148e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (146) ((1 180))) -1.813586922109707e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (147) ((1 180))) -1.6305767447865746e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (148) ((1 180))) -1.4571267217520587e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (149) ((1 180))) -1.2942542198392462e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (150) ((1 180))) -1.1424569194244598e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (151) ((1 180))) 0.0019282196424877589)
- (f2cl-lib:fset (f2cl-lib:fref alfa (152) ((1 180))) 0.0013559257630202223)
- (f2cl-lib:fset (f2cl-lib:fref alfa (153) ((1 180))) -7.178580904213031e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (154) ((1 180))) -0.0025808480257527035)
- (f2cl-lib:fset (f2cl-lib:fref alfa (155) ((1 180))) -0.003492711308261685)
- (f2cl-lib:fset (f2cl-lib:fref alfa (156) ((1 180))) -0.003469862993409606)
- (f2cl-lib:fset (f2cl-lib:fref alfa (157) ((1 180))) -0.002822852333513102)
- (f2cl-lib:fset (f2cl-lib:fref alfa (158) ((1 180))) -0.0018810307640489134)
- (f2cl-lib:fset (f2cl-lib:fref alfa (159) ((1 180))) -8.895317183839476e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (160) ((1 180))) 3.879121026310352e-6)
- (f2cl-lib:fset (f2cl-lib:fref alfa (161) ((1 180))) 7.286885401196914e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (162) ((1 180))) 0.0012656637305345775)
- (f2cl-lib:fset (f2cl-lib:fref alfa (163) ((1 180))) 0.0016251815837267443)
- (f2cl-lib:fset (f2cl-lib:fref alfa (164) ((1 180))) 0.0018320315321637317)
- (f2cl-lib:fset (f2cl-lib:fref alfa (165) ((1 180))) 0.001915883889905279)
- (f2cl-lib:fset (f2cl-lib:fref alfa (166) ((1 180))) 0.0019058884675554615)
- (f2cl-lib:fset (f2cl-lib:fref alfa (167) ((1 180))) 0.0018279898242182574)
- (f2cl-lib:fset (f2cl-lib:fref alfa (168) ((1 180))) 0.0017038950642112155)
- (f2cl-lib:fset (f2cl-lib:fref alfa (169) ((1 180))) 0.001550971271710977)
- (f2cl-lib:fset (f2cl-lib:fref alfa (170) ((1 180))) 0.0013826142185227616)
- (f2cl-lib:fset (f2cl-lib:fref alfa (171) ((1 180))) 0.0012088142423006478)
- (f2cl-lib:fset (f2cl-lib:fref alfa (172) ((1 180))) 0.0010367653263834496)
- (f2cl-lib:fset (f2cl-lib:fref alfa (173) ((1 180))) 8.714379180686191e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (174) ((1 180))) 7.16080155297701e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (175) ((1 180))) 5.726370025581293e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (176) ((1 180))) 4.420898194658023e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (177) ((1 180))) 3.247249485030906e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (178) ((1 180))) 2.2034204273024657e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (179) ((1 180))) 1.284128984013539e-4)
- (f2cl-lib:fset (f2cl-lib:fref alfa (180) ((1 180))) 4.8200592455209545e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (1) ((1 210))) 0.017998872141355333)
- (f2cl-lib:fset (f2cl-lib:fref beta (2) ((1 210))) 0.005599649110643881)
- (f2cl-lib:fset (f2cl-lib:fref beta (3) ((1 210))) 0.002885014022311328)
- (f2cl-lib:fset (f2cl-lib:fref beta (4) ((1 210))) 0.0018009660676105395)
- (f2cl-lib:fset (f2cl-lib:fref beta (5) ((1 210))) 0.001247531105891992)
- (f2cl-lib:fset (f2cl-lib:fref beta (6) ((1 210))) 9.228788765729384e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (7) ((1 210))) 7.144304217272874e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (8) ((1 210))) 5.71787281789705e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (9) ((1 210))) 4.6943100760648154e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (10) ((1 210))) 3.932328354629167e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (11) ((1 210))) 3.348188893182977e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (12) ((1 210))) 2.8895214849575157e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (13) ((1 210))) 2.5221161554957333e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (14) ((1 210))) 2.2228058079888333e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (15) ((1 210))) 1.9754183803306255e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (16) ((1 210))) 1.7683685501971802e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (17) ((1 210))) 1.593168996618211e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (18) ((1 210))) 1.44347930197334e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (19) ((1 210))) 1.314480681199654e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (20) ((1 210))) 1.2024544494930292e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (21) ((1 210))) 1.1044914450459939e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (22) ((1 210))) 1.0182877074056726e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (23) ((1 210))) 9.419982242042375e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (24) ((1 210))) 8.741305457538344e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (25) ((1 210))) 8.134662621628016e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (26) ((1 210))) 7.590022696462193e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (27) ((1 210))) 7.099063006341535e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (28) ((1 210))) 6.654828748424683e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (29) ((1 210))) 6.251469589692751e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (30) ((1 210))) 5.884033944262518e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (31) ((1 210))) -0.0014928295321342917)
- (f2cl-lib:fset (f2cl-lib:fref beta (32) ((1 210))) -8.782047095463893e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (33) ((1 210))) -5.029165495720346e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (34) ((1 210))) -2.948221385127461e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (35) ((1 210))) -1.7546399697078285e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (36) ((1 210))) -1.0400855046081643e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (37) ((1 210))) -5.96141953046458e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (38) ((1 210))) -3.120389290760983e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (39) ((1 210))) -1.2608973598023004e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (40) ((1 210))) -2.4289260857573036e-7)
- (f2cl-lib:fset (f2cl-lib:fref beta (41) ((1 210))) 8.059961654142735e-6)
- (f2cl-lib:fset (f2cl-lib:fref beta (42) ((1 210))) 1.3650700926214737e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (43) ((1 210))) 1.7396412547292628e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (44) ((1 210))) 1.9867297884213378e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (45) ((1 210))) 2.1446326379082264e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (46) ((1 210))) 2.2395465923245655e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (47) ((1 210))) 2.289677838147126e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (48) ((1 210))) 2.307853898111778e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (49) ((1 210))) 2.3032197608090915e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (50) ((1 210))) 2.2823607372034874e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (51) ((1 210))) 2.250058811052924e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (52) ((1 210))) 2.2098101536199147e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (53) ((1 210))) 2.1641842744810388e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (54) ((1 210))) 2.1150764925622084e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (55) ((1 210))) 2.0638874978217076e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (56) ((1 210))) 2.011652419970817e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (57) ((1 210))) 1.9591345014117925e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (58) ((1 210))) 1.9068936791043675e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (59) ((1 210))) 1.855337196416367e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (60) ((1 210))) 1.8047572225967418e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (61) ((1 210))) 5.522130767212928e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (62) ((1 210))) 4.4793258155238475e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (63) ((1 210))) 2.795206539920206e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (64) ((1 210))) 1.5246815619844659e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (65) ((1 210))) 6.932711056570436e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (66) ((1 210))) 1.7625868306999143e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (67) ((1 210))) -1.3574499634326914e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (68) ((1 210))) -3.1797241335042714e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (69) ((1 210))) -4.188618616966934e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (70) ((1 210))) -4.690048893791411e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (71) ((1 210))) -4.8766544741378737e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (72) ((1 210))) -4.8701003118673514e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (73) ((1 210))) -4.747556208900867e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (74) ((1 210))) -4.558130581386285e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (75) ((1 210))) -4.3330964451126613e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (76) ((1 210))) -4.0923019315775044e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (77) ((1 210))) -3.848226386032213e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (78) ((1 210))) -3.6085716753541053e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (79) ((1 210))) -3.3779330612336744e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (80) ((1 210))) -3.158885607721097e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (81) ((1 210))) -2.952695617508073e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (82) ((1 210))) -2.7597891482833575e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (83) ((1 210))) -2.580061746668837e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (84) ((1 210))) -2.413083567612802e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (85) ((1 210))) -2.2582350951834607e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (86) ((1 210))) -2.1147965676891295e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (87) ((1 210))) -1.982006388852949e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (88) ((1 210))) -1.8590987080106508e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (89) ((1 210))) -1.7453269984421022e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (90) ((1 210))) -1.6399782385449801e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (91) ((1 210))) -4.746177965599597e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (92) ((1 210))) -4.778645671473215e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (93) ((1 210))) -3.2039022806703765e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (94) ((1 210))) -1.6110501611996228e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (95) ((1 210))) -4.2577810128543525e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (96) ((1 210))) 3.445712942949675e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (97) ((1 210))) 7.97092684075675e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (98) ((1 210))) 1.031382367082722e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (99) ((1 210))) 1.1246677526220417e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (100) ((1 210))) 1.1310364210848138e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (101) ((1 210))) 1.0865163484877427e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (102) ((1 210))) 1.0143795159766197e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (103) ((1 210))) 9.292983965933639e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (104) ((1 210))) 8.402931330160901e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (105) ((1 210))) 7.527279913491342e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (106) ((1 210))) 6.696325219757311e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (107) ((1 210))) 5.925645473231947e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (108) ((1 210))) 5.221693088269757e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (109) ((1 210))) 4.585394851653607e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (110) ((1 210))) 4.014455138914868e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (111) ((1 210))) 3.5048173003132804e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (112) ((1 210))) 3.0515799503434665e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (113) ((1 210))) 2.6495611995051604e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (114) ((1 210))) 2.293636336909982e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (115) ((1 210))) 1.9789305666402165e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (116) ((1 210))) 1.7009198463641262e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (117) ((1 210))) 1.4554742826152403e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (118) ((1 210))) 1.2388664099587843e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (119) ((1 210))) 1.0477587607658323e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (120) ((1 210))) 8.791799549784793e-6)
- (f2cl-lib:fset (f2cl-lib:fref beta (121) ((1 210))) 7.364658105725783e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (122) ((1 210))) 8.72790805146194e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (123) ((1 210))) 6.22614862573135e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (124) ((1 210))) 2.859981541943042e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (125) ((1 210))) 3.8473767287936605e-6)
- (f2cl-lib:fset (f2cl-lib:fref beta (126) ((1 210))) -1.8790600363697157e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (127) ((1 210))) -2.9760364659455457e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (128) ((1 210))) -3.459981268326564e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (129) ((1 210))) -3.533824709160377e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (130) ((1 210))) -3.357156357750488e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (131) ((1 210))) -3.043211247890398e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (132) ((1 210))) -2.667227230476128e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (133) ((1 210))) -2.2765421412281958e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (134) ((1 210))) -1.8992261185456238e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (135) ((1 210))) -1.5505891859909388e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (136) ((1 210))) -1.2377824076187363e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (137) ((1 210))) -9.629261477176442e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (138) ((1 210))) -7.251783277144254e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (139) ((1 210))) -5.220700288956338e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (140) ((1 210))) -3.5034775051190054e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (141) ((1 210))) -2.0648976103555175e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (142) ((1 210))) -8.70106096849767e-6)
- (f2cl-lib:fset (f2cl-lib:fref beta (143) ((1 210))) 1.1369868667510027e-6)
- (f2cl-lib:fset (f2cl-lib:fref beta (144) ((1 210))) 9.164264741227788e-6)
- (f2cl-lib:fset (f2cl-lib:fref beta (145) ((1 210))) 1.5647778542887264e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (146) ((1 210))) 2.0822362948246687e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (147) ((1 210))) 2.4892338100459516e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (148) ((1 210))) 2.8034050957414636e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (149) ((1 210))) 3.0398777462986193e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (150) ((1 210))) 3.2115673140670065e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (151) ((1 210))) -0.0018018219196388571)
- (f2cl-lib:fset (f2cl-lib:fref beta (152) ((1 210))) -0.0024340296293804253)
- (f2cl-lib:fset (f2cl-lib:fref beta (153) ((1 210))) -0.0018342266354985683)
- (f2cl-lib:fset (f2cl-lib:fref beta (154) ((1 210))) -7.622045963540097e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (155) ((1 210))) 2.3907947525692724e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (156) ((1 210))) 9.49266117176881e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (157) ((1 210))) 0.0013446744970154036)
- (f2cl-lib:fset (f2cl-lib:fref beta (158) ((1 210))) 0.0014845749525944918)
- (f2cl-lib:fset (f2cl-lib:fref beta (159) ((1 210))) 0.001447323398306176)
- (f2cl-lib:fset (f2cl-lib:fref beta (160) ((1 210))) 0.0013026826128565718)
- (f2cl-lib:fset (f2cl-lib:fref beta (161) ((1 210))) 0.0011035159737564268)
- (f2cl-lib:fset (f2cl-lib:fref beta (162) ((1 210))) 8.860474404197918e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (163) ((1 210))) 6.730732081656655e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (164) ((1 210))) 4.776038728565824e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (165) ((1 210))) 3.0599192635878936e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (166) ((1 210))) 1.6031569459472167e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (167) ((1 210))) 4.0074955527061334e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (168) ((1 210))) -5.666074616352517e-5)
- (f2cl-lib:fset (f2cl-lib:fref beta (169) ((1 210))) -1.3250618677298265e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (170) ((1 210))) -1.9029618798961406e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (171) ((1 210))) -2.328114503769374e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (172) ((1 210))) -2.626288114646688e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (173) ((1 210))) -2.8205046986759863e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (174) ((1 210))) -2.9308156319286116e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (175) ((1 210))) -2.974359621763166e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (176) ((1 210))) -2.9655733423934805e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (177) ((1 210))) -2.9164736331209085e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (178) ((1 210))) -2.836962038377342e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (179) ((1 210))) -2.7351231709567336e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (180) ((1 210))) -2.617501558067686e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (181) ((1 210))) 0.006385858912120509)
- (f2cl-lib:fset (f2cl-lib:fref beta (182) ((1 210))) 0.009623742158063778)
- (f2cl-lib:fset (f2cl-lib:fref beta (183) ((1 210))) 0.007618780612070011)
- (f2cl-lib:fset (f2cl-lib:fref beta (184) ((1 210))) 0.002832190555456281)
- (f2cl-lib:fset (f2cl-lib:fref beta (185) ((1 210))) -0.002098413520127201)
- (f2cl-lib:fset (f2cl-lib:fref beta (186) ((1 210))) -0.005738267642166265)
- (f2cl-lib:fset (f2cl-lib:fref beta (187) ((1 210))) -0.0077080424449541465)
- (f2cl-lib:fset (f2cl-lib:fref beta (188) ((1 210))) -0.008210116922648444)
- (f2cl-lib:fset (f2cl-lib:fref beta (189) ((1 210))) -0.007658245203469054)
- (f2cl-lib:fset (f2cl-lib:fref beta (190) ((1 210))) -0.006472097293910452)
- (f2cl-lib:fset (f2cl-lib:fref beta (191) ((1 210))) -0.004991324120049665)
- (f2cl-lib:fset (f2cl-lib:fref beta (192) ((1 210))) -0.003456122897131333)
- (f2cl-lib:fset (f2cl-lib:fref beta (193) ((1 210))) -0.002017855800141708)
- (f2cl-lib:fset (f2cl-lib:fref beta (194) ((1 210))) -7.594306867819615e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (195) ((1 210))) 2.8417363152385916e-4)
- (f2cl-lib:fset (f2cl-lib:fref beta (196) ((1 210))) 0.0011089166758633739)
- (f2cl-lib:fset (f2cl-lib:fref beta (197) ((1 210))) 0.0017290149387272878)
- (f2cl-lib:fset (f2cl-lib:fref beta (198) ((1 210))) 0.0021681259080268468)
- (f2cl-lib:fset (f2cl-lib:fref beta (199) ((1 210))) 0.0024535771049453976)
- (f2cl-lib:fset (f2cl-lib:fref beta (200) ((1 210))) 0.0026128182105833483)
- (f2cl-lib:fset (f2cl-lib:fref beta (201) ((1 210))) 0.002671410396562769)
- (f2cl-lib:fset (f2cl-lib:fref beta (202) ((1 210))) 0.0026520307339598045)
- (f2cl-lib:fset (f2cl-lib:fref beta (203) ((1 210))) 0.002574116528772873)
- (f2cl-lib:fset (f2cl-lib:fref beta (204) ((1 210))) 0.0024538912623609443)
- (f2cl-lib:fset (f2cl-lib:fref beta (205) ((1 210))) 0.002304600580717955)
- (f2cl-lib:fset (f2cl-lib:fref beta (206) ((1 210))) 0.0021368483768671267)
- (f2cl-lib:fset (f2cl-lib:fref beta (207) ((1 210))) 0.001958965284788709)
- (f2cl-lib:fset (f2cl-lib:fref beta (208) ((1 210))) 0.0017773700867945441)
- (f2cl-lib:fset (f2cl-lib:fref beta (209) ((1 210))) 0.0015969028076583908)
- (f2cl-lib:fset (f2cl-lib:fref beta (210) ((1 210))) 0.0014211197566443854)
- (f2cl-lib:fset (f2cl-lib:fref gama (1) ((1 30))) 0.6299605249474367)
- (f2cl-lib:fset (f2cl-lib:fref gama (2) ((1 30))) 0.25198420997897464)
- (f2cl-lib:fset (f2cl-lib:fref gama (3) ((1 30))) 0.15479030041565586)
- (f2cl-lib:fset (f2cl-lib:fref gama (4) ((1 30))) 0.11071306241615902)
- (f2cl-lib:fset (f2cl-lib:fref gama (5) ((1 30))) 0.08573093955273949)
- (f2cl-lib:fset (f2cl-lib:fref gama (6) ((1 30))) 0.06971613169586843)
- (f2cl-lib:fset (f2cl-lib:fref gama (7) ((1 30))) 0.05860856718937136)
- (f2cl-lib:fset (f2cl-lib:fref gama (8) ((1 30))) 0.050469887353631074)
- (f2cl-lib:fset (f2cl-lib:fref gama (9) ((1 30))) 0.04426005806891548)
- (f2cl-lib:fset (f2cl-lib:fref gama (10) ((1 30))) 0.039372066154350994)
- (f2cl-lib:fset (f2cl-lib:fref gama (11) ((1 30))) 0.03542831959244554)
- (f2cl-lib:fset (f2cl-lib:fref gama (12) ((1 30))) 0.032181885750209825)
- (f2cl-lib:fset (f2cl-lib:fref gama (13) ((1 30))) 0.029464624079115768)
- (f2cl-lib:fset (f2cl-lib:fref gama (14) ((1 30))) 0.027158167711293448)
- (f2cl-lib:fset (f2cl-lib:fref gama (15) ((1 30))) 0.025176827297386177)
- (f2cl-lib:fset (f2cl-lib:fref gama (16) ((1 30))) 0.02345707553060789)
- (f2cl-lib:fset (f2cl-lib:fref gama (17) ((1 30))) 0.021950839013490722)
- (f2cl-lib:fset (f2cl-lib:fref gama (18) ((1 30))) 0.020621082823564625)
- (f2cl-lib:fset (f2cl-lib:fref gama (19) ((1 30))) 0.019438824089788084)
- (f2cl-lib:fset (f2cl-lib:fref gama (20) ((1 30))) 0.018381063380068317)
- (f2cl-lib:fset (f2cl-lib:fref gama (21) ((1 30))) 0.017429321323196318)
- (f2cl-lib:fset (f2cl-lib:fref gama (22) ((1 30))) 0.016568583778661234)
- (f2cl-lib:fset (f2cl-lib:fref gama (23) ((1 30))) 0.015786528598791844)
- (f2cl-lib:fset (f2cl-lib:fref gama (24) ((1 30))) 0.01507295014940956)
- (f2cl-lib:fset (f2cl-lib:fref gama (25) ((1 30))) 0.014419325083995465)
- (f2cl-lib:fset (f2cl-lib:fref gama (26) ((1 30))) 0.013818480573534178)
- (f2cl-lib:fset (f2cl-lib:fref gama (27) ((1 30))) 0.013264337899427657)
- (f2cl-lib:fset (f2cl-lib:fref gama (28) ((1 30))) 0.012751712197049864)
- (f2cl-lib:fset (f2cl-lib:fref gama (29) ((1 30))) 0.012276154531876278)
- (f2cl-lib:fset (f2cl-lib:fref gama (30) ((1 30))) 0.01183382623984824)
- (defun zunhj
- (zr zi fnu ipmtr tol phir phii argr argi zeta1r zeta1i zeta2r zeta2i
- asumr asumi bsumr bsumi)
- (declare (type f2cl-lib:integer4 ipmtr)
- (type double-float zr zi fnu tol phir phii argr argi zeta1r zeta1i
- zeta2r zeta2i asumr asumi bsumr bsumi))
- (prog ((ap (make-array 30 :element-type 'double-float))
- (pr (make-array 30 :element-type 'double-float))
- (pi_ (make-array 30 :element-type 'double-float))
- (upr (make-array 14 :element-type 'double-float))
- (upi (make-array 14 :element-type 'double-float))
- (crr (make-array 14 :element-type 'double-float))
- (cri (make-array 14 :element-type 'double-float))
- (drr (make-array 14 :element-type 'double-float))
- (dri (make-array 14 :element-type 'double-float)) (ias 0) (ibs 0)
- (is 0) (j 0) (jr 0) (ju 0) (k 0) (kmax 0) (kp1 0) (ks 0) (l 0)
- (lr 0) (lrp1 0) (l1 0) (l2 0) (m 0) (idum 0) (ang 0.0) (atol 0.0)
- (aw2 0.0) (azth 0.0) (btol 0.0) (fn13 0.0) (fn23 0.0) (pp 0.0)
- (przthi 0.0) (przthr 0.0) (ptfni 0.0) (ptfnr 0.0) (raw 0.0)
- (raw2 0.0) (razth 0.0) (rfnu 0.0) (rfnu2 0.0) (rfn13 0.0)
- (rtzti 0.0) (rtztr 0.0) (rzthi 0.0) (rzthr 0.0) (sti 0.0) (str 0.0)
- (sumai 0.0) (sumar 0.0) (sumbi 0.0) (sumbr 0.0) (test 0.0)
- (tfni 0.0) (tfnr 0.0) (tzai 0.0) (tzar 0.0) (t2i 0.0) (t2r 0.0)
- (wi 0.0) (wr 0.0) (w2i 0.0) (w2r 0.0) (zai 0.0) (zar 0.0) (zbi 0.0)
- (zbr 0.0) (zci 0.0) (zcr 0.0) (zetai 0.0) (zetar 0.0) (zthi 0.0)
- (zthr 0.0) (ac 0.0))
- (declare (type (simple-array double-float (14)) upr upi drr dri crr cri)
- (type (simple-array double-float (30)) pr pi_ ap)
- (type double-float ac zthr zthi zetar zetai zcr zci zbr zbi zar
- zai w2r w2i wr wi t2r t2i tzar tzai tfnr tfni test sumbr sumbi
- sumar sumai str sti rzthr rzthi rtztr rtzti rfn13 rfnu2 rfnu
- razth raw2 raw ptfnr ptfni przthr przthi pp fn23 fn13 btol azth
- aw2 atol ang)
- (type f2cl-lib:integer4 idum m l2 l1 lrp1 lr l ks kp1 kmax k ju
- jr j is ibs ias))
- (setf rfnu (/ 1.0 fnu))
- (setf test (* (f2cl-lib:d1mach 1) 1000.0))
- (setf ac (* fnu test))
- (if (or (> (abs zr) ac) (> (abs zi) ac)) (go label15))
- (setf zeta1r (+ (* 2.0 (abs (f2cl-lib:flog test))) fnu))
- (setf zeta1i 0.0)
- (setf zeta2r fnu)
- (setf zeta2i 0.0)
- (setf phir 1.0)
- (setf phii 0.0)
- (setf argr 1.0)
- (setf argi 0.0)
- (go end_label)
- label15
- (setf zbr (* zr rfnu))
- (setf zbi (* zi rfnu))
- (setf rfnu2 (* rfnu rfnu))
- (setf fn13 (expt fnu ex1))
- (setf fn23 (* fn13 fn13))
- (setf rfn13 (/ 1.0 fn13))
- (setf w2r (+ (- coner (* zbr zbr)) (* zbi zbi)))
- (setf w2i (- conei (* zbr zbi) (* zbr zbi)))
- (setf aw2 (zabs w2r w2i))
- (if (> aw2 0.25) (go label130))
- (setf k 1)
- (f2cl-lib:fset (f2cl-lib:fref pr (1) ((1 30))) coner)
- (f2cl-lib:fset (f2cl-lib:fref pi_ (1) ((1 30))) conei)
- (setf sumar (f2cl-lib:fref gama (1) ((1 30))))
- (setf sumai zeroi)
- (f2cl-lib:fset (f2cl-lib:fref ap (1) ((1 30))) 1.0)
- (if (< aw2 tol) (go label20))
- (f2cl-lib:fdo (k 2 (f2cl-lib:int-add k 1))
- ((> k 30) nil)
- (tagbody
- (f2cl-lib:fset (f2cl-lib:fref pr (k) ((1 30)))
- (-
- (*
- (f2cl-lib:fref pr ((f2cl-lib:int-sub k 1)) ((1 30)))
- w2r)
- (*
- (f2cl-lib:fref pi_
- ((f2cl-lib:int-sub k 1))
- ((1 30)))
- w2i)))
- (f2cl-lib:fset (f2cl-lib:fref pi_ (k) ((1 30)))
- (+
- (*
- (f2cl-lib:fref pr ((f2cl-lib:int-sub k 1)) ((1 30)))
- w2i)
- (*
- (f2cl-lib:fref pi_
- ((f2cl-lib:int-sub k 1))
- ((1 30)))
- w2r)))
- (setf sumar
- (+ sumar
- (* (f2cl-lib:fref pr (k) ((1 30)))
- (f2cl-lib:fref gama (k) ((1 30))))))
- (setf sumai
- (+ sumai
- (* (f2cl-lib:fref pi_ (k) ((1 30)))
- (f2cl-lib:fref gama (k) ((1 30))))))
- (f2cl-lib:fset (f2cl-lib:fref ap (k) ((1 30)))
- (*
- (f2cl-lib:fref ap ((f2cl-lib:int-sub k 1)) ((1 30)))
- aw2))
- (if (< (f2cl-lib:fref ap (k) ((1 30))) tol) (go label20))
- label10))
- (setf k 30)
- label20
- (setf kmax k)
- (setf zetar (- (* w2r sumar) (* w2i sumai)))
- (setf zetai (+ (* w2r sumai) (* w2i sumar)))
- (setf argr (* zetar fn23))
- (setf argi (* zetai fn23))
- (multiple-value-bind
- (var-0 var-1 var-2 var-3)
- (zsqrt sumar sumai zar zai)
- (declare (ignore var-0 var-1))
- (setf zar var-2)
- (setf zai var-3))
- (multiple-value-bind
- (var-0 var-1 var-2 var-3)
- (zsqrt w2r w2i str sti)
- (declare (ignore var-0 var-1))
- (setf str var-2)
- (setf sti var-3))
- (setf zeta2r (* str fnu))
- (setf zeta2i (* sti fnu))
- (setf str (+ coner (* ex2 (- (* zetar zar) (* zetai zai)))))
- (setf sti (+ conei (* ex2 (+ (* zetar zai) (* zetai zar)))))
- (setf zeta1r (- (* str zeta2r) (* sti zeta2i)))
- (setf zeta1i (+ (* str zeta2i) (* sti zeta2r)))
- (setf zar (+ zar zar))
- (setf zai (+ zai zai))
- (multiple-value-bind
- (var-0 var-1 var-2 var-3)
- (zsqrt zar zai str sti)
- (declare (ignore var-0 var-1))
- (setf str var-2)
- (setf sti var-3))
- (setf phir (* str rfn13))
- (setf phii (* sti rfn13))
- (if (= ipmtr 1) (go label120))
- (setf sumbr zeror)
- (setf sumbi zeroi)
- (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
- ((> k kmax) nil)
- (tagbody
- (setf sumbr
- (+ sumbr
- (* (f2cl-lib:fref pr (k) ((1 30)))
- (f2cl-lib:fref beta (k) ((1 210))))))
- (setf sumbi
- (+ sumbi
- (* (f2cl-lib:fref pi_ (k) ((1 30)))
- (f2cl-lib:fref beta (k) ((1 210))))))
- label30))
- (setf asumr zeror)
- (setf asumi zeroi)
- (setf bsumr sumbr)
- (setf bsumi sumbi)
- (setf l1 0)
- (setf l2 30)
- (setf btol (* tol (+ (abs bsumr) (abs bsumi))))
- (setf atol tol)
- (setf pp 1.0)
- (setf ias 0)
- (setf ibs 0)
- (if (< rfnu2 tol) (go label110))
- (f2cl-lib:fdo (is 2 (f2cl-lib:int-add is 1))
- ((> is 7) nil)
- (tagbody
- (setf atol (/ atol rfnu2))
- (setf pp (* pp rfnu2))
- (if (= ias 1) (go label60))
- (setf sumar zeror)
- (setf sumai zeroi)
- (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
- ((> k kmax) nil)
- (tagbody
- (setf m (f2cl-lib:int-add l1 k))
- (setf sumar
- (+ sumar
- (* (f2cl-lib:fref pr (k) ((1 30)))
- (f2cl-lib:fref alfa (m) ((1 180))))))
- (setf sumai
- (+ sumai
- (* (f2cl-lib:fref pi_ (k) ((1 30)))
- (f2cl-lib:fref alfa (m) ((1 180))))))
- (if (< (f2cl-lib:fref ap (k) ((1 30))) atol) (go label50))
- label40))
- label50
- (setf asumr (+ asumr (* sumar pp)))
- (setf asumi (+ asumi (* sumai pp)))
- (if (< pp tol) (setf ias 1))
- label60
- (if (= ibs 1) (go label90))
- (setf sumbr zeror)
- (setf sumbi zeroi)
- (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
- ((> k kmax) nil)
- (tagbody
- (setf m (f2cl-lib:int-add l2 k))
- (setf sumbr
- (+ sumbr
- (* (f2cl-lib:fref pr (k) ((1 30)))
- (f2cl-lib:fref beta (m) ((1 210))))))
- (setf sumbi
- (+ sumbi
- (* (f2cl-lib:fref pi_ (k) ((1 30)))
- (f2cl-lib:fref beta (m) ((1 210))))))
- (if (< (f2cl-lib:fref ap (k) ((1 30))) atol) (go label80))
- label70))
- label80
- (setf bsumr (+ bsumr (* sumbr pp)))
- (setf bsumi (+ bsumi (* sumbi pp)))
- (if (< pp btol) (setf ibs 1))
- label90
- (if (and (= ias 1) (= ibs 1)) (go label110))
- (setf l1 (f2cl-lib:int-add l1 30))
- (setf l2 (f2cl-lib:int-add l2 30))
- label100))
- label110
- (setf asumr (+ asumr coner))
- (setf pp (* rfnu rfn13))
- (setf bsumr (* bsumr pp))
- (setf bsumi (* bsumi pp))
- label120
- (go end_label)
- label130
- (multiple-value-bind
- (var-0 var-1 var-2 var-3)
- (zsqrt w2r w2i wr wi)
- (declare (ignore var-0 var-1))
- (setf wr var-2)
- (setf wi var-3))
- (if (< wr 0.0) (setf wr 0.0))
- (if (< wi 0.0) (setf wi 0.0))
- (setf str (+ coner wr))
- (setf sti wi)
- (multiple-value-bind
- (var-0 var-1 var-2 var-3 var-4 var-5)
- (zdiv str sti zbr zbi zar zai)
- (declare (ignore var-0 var-1 var-2 var-3))
- (setf zar var-4)
- (setf zai var-5))
- (multiple-value-bind
- (var-0 var-1 var-2 var-3 var-4)
- (zlog zar zai zcr zci idum)
- (declare (ignore var-0 var-1))
- (setf zcr var-2)
- (setf zci var-3)
- (setf idum var-4))
- (if (< zci 0.0) (setf zci 0.0))
- (if (> zci hpi) (setf zci hpi))
- (if (< zcr 0.0) (setf zcr 0.0))
- (setf zthr (* (- zcr wr) 1.5))
- (setf zthi (* (- zci wi) 1.5))
- (setf zeta1r (* zcr fnu))
- (setf zeta1i (* zci fnu))
- (setf zeta2r (* wr fnu))
- (setf zeta2i (* wi fnu))
- (setf azth (zabs zthr zthi))
- (setf ang thpi)
- (if (and (>= zthr 0.0) (< zthi 0.0)) (go label140))
- (setf ang hpi)
- (if (= zthr 0.0) (go label140))
- (setf ang (f2cl-lib:datan (/ zthi zthr)))
- (if (< zthr 0.0) (setf ang (+ ang gpi)))
- label140
- (setf pp (expt azth ex2))
- (setf ang (* ang ex2))
- (setf zetar (* pp (cos ang)))
- (setf zetai (* pp (sin ang)))
- (if (< zetai 0.0) (setf zetai 0.0))
- (setf argr (* zetar fn23))
- (setf argi (* zetai fn23))
- (multiple-value-bind
- (var-0 var-1 var-2 var-3 var-4 var-5)
- (zdiv zthr zthi zetar zetai rtztr rtzti)
- (declare (ignore var-0 var-1 var-2 var-3))
- (setf rtztr var-4)
- (setf rtzti var-5))
- (multiple-value-bind
- (var-0 var-1 var-2 var-3 var-4 var-5)
- (zdiv rtztr rtzti wr wi zar zai)
- (declare (ignore var-0 var-1 var-2 var-3))
- (setf zar var-4)
- (setf zai var-5))
- (setf tzar (+ zar zar))
- (setf tzai (+ zai zai))
- (multiple-value-bind
- (var-0 var-1 var-2 var-3)
- (zsqrt tzar tzai str sti)
- (declare (ignore var-0 var-1))
- (setf str var-2)
- (setf sti var-3))
- (setf phir (* str rfn13))
- (setf phii (* sti rfn13))
- (if (= ipmtr 1) (go label120))
- (setf raw (/ 1.0 (f2cl-lib:fsqrt aw2)))
- (setf str (* wr raw))
- (setf sti (* (- wi) raw))
- (setf tfnr (* str rfnu raw))
- (setf tfni (* sti rfnu raw))
- (setf razth (/ 1.0 azth))
- (setf str (* zthr razth))
- (setf sti (* (- zthi) razth))
- (setf rzthr (* str razth rfnu))
- (setf rzthi (* sti razth rfnu))
- (setf zcr (* rzthr (f2cl-lib:fref ar (2) ((1 14)))))
- (setf zci (* rzthi (f2cl-lib:fref ar (2) ((1 14)))))
- (setf raw2 (/ 1.0 aw2))
- (setf str (* w2r raw2))
- (setf sti (* (- w2i) raw2))
- (setf t2r (* str raw2))
- (setf t2i (* sti raw2))
- (setf str
- (+ (* t2r (f2cl-lib:fref c (2) ((1 105))))
- (f2cl-lib:fref c (3) ((1 105)))))
- (setf sti (* t2i (f2cl-lib:fref c (2) ((1 105)))))
- (f2cl-lib:fset (f2cl-lib:fref upr (2) ((1 14)))
- (- (* str tfnr) (* sti tfni)))
- (f2cl-lib:fset (f2cl-lib:fref upi (2) ((1 14)))
- (+ (* str tfni) (* sti tfnr)))
- (setf bsumr (+ (f2cl-lib:fref upr (2) ((1 14))) zcr))
- (setf bsumi (+ (f2cl-lib:fref upi (2) ((1 14))) zci))
- (setf asumr zeror)
- (setf asumi zeroi)
- (if (< rfnu tol) (go label220))
- (setf przthr rzthr)
- (setf przthi rzthi)
- (setf ptfnr tfnr)
- (setf ptfni tfni)
- (f2cl-lib:fset (f2cl-lib:fref upr (1) ((1 14))) coner)
- (f2cl-lib:fset (f2cl-lib:fref upi (1) ((1 14))) conei)
- (setf pp 1.0)
- (setf btol (* tol (+ (abs bsumr) (abs bsumi))))
- (setf ks 0)
- (setf kp1 2)
- (setf l 3)
- (setf ias 0)
- (setf ibs 0)
- (f2cl-lib:fdo (lr 2 (f2cl-lib:int-add lr 2))
- ((> lr 12) nil)
- (tagbody
- (setf lrp1 (f2cl-lib:int-add lr 1))
- (f2cl-lib:fdo (k lr (f2cl-lib:int-add k 1))
- ((> k lrp1) nil)
- (tagbody
- (setf ks (f2cl-lib:int-add ks 1))
- (setf kp1 (f2cl-lib:int-add kp1 1))
- (setf l (f2cl-lib:int-add l 1))
- (setf zar (f2cl-lib:fref c (l) ((1 105))))
- (setf zai zeroi)
- (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
- ((> j kp1) nil)
- (tagbody
- (setf l (f2cl-lib:int-add l 1))
- (setf str
- (+ (- (* zar t2r) (* t2i zai))
- (f2cl-lib:fref c (l) ((1 105)))))
- (setf zai (+ (* zar t2i) (* zai t2r)))
- (setf zar str)
- label150))
- (setf str (- (* ptfnr tfnr) (* ptfni tfni)))
- (setf ptfni (+ (* ptfnr tfni) (* ptfni tfnr)))
- (setf ptfnr str)
- (f2cl-lib:fset (f2cl-lib:fref upr (kp1) ((1 14)))
- (- (* ptfnr zar) (* ptfni zai)))
- (f2cl-lib:fset (f2cl-lib:fref upi (kp1) ((1 14)))
- (+ (* ptfni zar) (* ptfnr zai)))
- (f2cl-lib:fset (f2cl-lib:fref crr (ks) ((1 14)))
- (* przthr
- (f2cl-lib:fref br
- ((f2cl-lib:int-add ks 1))
- ((1 14)))))
- (f2cl-lib:fset (f2cl-lib:fref cri (ks) ((1 14)))
- (* przthi
- (f2cl-lib:fref br
- ((f2cl-lib:int-add ks 1))
- ((1 14)))))
- (setf str (- (* przthr rzthr) (* przthi rzthi)))
- (setf przthi (+ (* przthr rzthi) (* przthi rzthr)))
- (setf przthr str)
- (f2cl-lib:fset (f2cl-lib:fref drr (ks) ((1 14)))
- (* przthr
- (f2cl-lib:fref ar
- ((f2cl-lib:int-add ks 2))
- ((1 14)))))
- (f2cl-lib:fset (f2cl-lib:fref dri (ks) ((1 14)))
- (* przthi
- (f2cl-lib:fref ar
- ((f2cl-lib:int-add ks 2))
- ((1 14)))))
- label160))
- (setf pp (* pp rfnu2))
- (if (= ias 1) (go label180))
- (setf sumar (f2cl-lib:fref upr (lrp1) ((1 14))))
- (setf sumai (f2cl-lib:fref upi (lrp1) ((1 14))))
- (setf ju lrp1)
- (f2cl-lib:fdo (jr 1 (f2cl-lib:int-add jr 1))
- ((> jr lr) nil)
- (tagbody
- (setf ju (f2cl-lib:int-sub ju 1))
- (setf sumar
- (-
- (+ sumar
- (* (f2cl-lib:fref crr (jr) ((1 14)))
- (f2cl-lib:fref upr (ju) ((1 14)))))
- (* (f2cl-lib:fref cri (jr) ((1 14)))
- (f2cl-lib:fref upi (ju) ((1 14))))))
- (setf sumai
- (+ sumai
- (* (f2cl-lib:fref crr (jr) ((1 14)))
- (f2cl-lib:fref upi (ju) ((1 14))))
- (* (f2cl-lib:fref cri (jr) ((1 14)))
- (f2cl-lib:fref upr (ju) ((1 14))))))
- label170))
- (setf asumr (+ asumr sumar))
- (setf asumi (+ asumi sumai))
- (setf test (coerce (+ (abs sumar) (abs sumai)) 'double-float))
- (if (and (< pp tol) (< test tol)) (setf ias 1))
- label180
- (if (= ibs 1) (go label200))
- (setf sumbr
- (-
- (+ (f2cl-lib:fref upr ((f2cl-lib:int-add lr 2)) ((1 14)))
- (* (f2cl-lib:fref upr (lrp1) ((1 14))) zcr))
- (* (f2cl-lib:fref upi (lrp1) ((1 14))) zci)))
- (setf sumbi
- (+ (f2cl-lib:fref upi ((f2cl-lib:int-add lr 2)) ((1 14)))
- (* (f2cl-lib:fref upr (lrp1) ((1 14))) zci)
- (* (f2cl-lib:fref upi (lrp1) ((1 14))) zcr)))
- (setf ju lrp1)
- (f2cl-lib:fdo (jr 1 (f2cl-lib:int-add jr 1))
- ((> jr lr) nil)
- (tagbody
- (setf ju (f2cl-lib:int-sub ju 1))
- (setf sumbr
- (-
- (+ sumbr
- (* (f2cl-lib:fref drr (jr) ((1 14)))
- (f2cl-lib:fref upr (ju) ((1 14)))))
- (* (f2cl-lib:fref dri (jr) ((1 14)))
- (f2cl-lib:fref upi (ju) ((1 14))))))
- (setf sumbi
- (+ sumbi
- (* (f2cl-lib:fref drr (jr) ((1 14)))
- (f2cl-lib:fref upi (ju) ((1 14))))
- (* (f2cl-lib:fref dri (jr) ((1 14)))
- (f2cl-lib:fref upr (ju) ((1 14))))))
- label190))
- (setf bsumr (+ bsumr sumbr))
- (setf bsumi (+ bsumi sumbi))
- (setf test (coerce (+ (abs sumbr) (abs sumbi)) 'double-float))
- (if (and (< pp btol) (< test btol)) (setf ibs 1))
- label200
- (if (and (= ias 1) (= ibs 1)) (go label220))
- label210))
- label220
- (setf asumr (+ asumr coner))
- (setf str (* (- bsumr) rfn13))
- (setf sti (* (- bsumi) rfn13))
- (multiple-value-bind
- (var-0 var-1 var-2 var-3 var-4 var-5)
- (zdiv str sti rtztr rtzti bsumr bsumi)
- (declare (ignore var-0 var-1 var-2 var-3))
- (setf bsumr var-4)
- (setf bsumi var-5))
- (go label120)
- end_label
- (return
- (values nil
- nil
- nil
- nil
- nil
- phir
- phii
- argr
- argi
- zeta1r
- zeta1i
- zeta2r
- zeta2i
- asumr
- asumi
- bsumr
- bsumi)))))
-
-